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The splitting of neutrinos (i^ — > wv) , a viable reaction for superluminal neutrinos, is shown to have 
phenomenologically relevant consequences if one accepts the recent OPERA results for the velocity of 
neutrinos. Neutrino splitting can be used to put strong constraints on the energy dependence of the 
velocity of propagation of neutrinos in a general analysis of modifications of relativistic kinematics 
and to propose observable effects due to the departures from special relativity in neutrino physics. 
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INTRODUCTION 

Neutrinos are a special window to physics beyond the 
Standard Model [H . They are the only particles which do 
not experience any long range interaction, which makes 
them particularly sensitive to possible new physics at 
the TeV scale or beyond. One exarnple is the possible 
existence of large extra-dimensions [2| which could af- 
fect neutrinos but not the remaining elementary parti- 
cles. Another example is the related property of neu- 
trinos of being able to propagate freely over very large 
distances. This provides an amplification mechanism for 
effects coming from new physics. 

The special role of neutrinos has been considered also 
in connection with possible departures from special rela- 
tivity (SR) and CPT invariance This can be under- 
stood if neutrinos have a special ability to explore the 
quantum nature of space-time, which is the more plau- 
sible origin for a departure from SR In this way, 
corrections to the relativistic kinematics for neutrinos 
have been considered previously as a possible alterna- 
tive (or coniplement to masses) explanation of neutrino 
oscillations [5[ . Very high- neutrino physics has also been 
suggested as a tool to examine quantum-gravity-induced 
corrections to the relativistic dispersion relation Q . 

Recently the OPERA neutrino experiment has re- 
ported the observation of a deviation in the time 
of flight of neutrinos with respect to SR correspond- 
ing to a relative difference of the muon neutrino ve- 
locity with respect to the speed of light of^ f — 1 ~ 
(2.48 ± 0.28 (stat.) ± 0.30 (sys.)) x 10"^ This rcsuh, 
together with similar although less statistically signifi- 
cant indications in the same direction from another long 
base line experiment (MINOS) Q, comes out as a major 
surprise since the corrections to SR it implies are much 
larger than what one would expect if they came from the 
natural scale of quantum gravity Q. 

What is more, this result is in contrast with the strin- 
gent limit oi V — 1 < 2 x 10~^ set by the observation of 
(anti) neutrinos emitted by the SN1987A supernova [loj . 
Although we are dealing here with a different type of neu- 



^ We use units where the speed of light c is equal to one. 



trinos (muon neutrinos in the accelerator experiments, 
versus electron neutrinos in the supernova case), fiavor- 
dependcnt Lorentz violations are severely constrained 
by neutrino oscillations However, neutrinos from 

SN1987A arc in the 10 MeV range, three orders of mag- 
nitude below the energy involved in the OPERA/MINOS 
experiments, which suggests the possibility that the ap- 
parent confiict between these observations might be due 
to an energy dependence of the deviations from SR, which 
is also natural from the point of view of the emergence 
of new physics controlled by a high-energy scale. All this 
provides a further motivation to consider the phcnomeno- 
logical consistency of superluminal neutrinos from a dif- 
ferent perspective. 

Among the observations which are sensitive to depar- 
tures from SR a special role is played by those related 
to changes of thresholds and those related to the con- 
version of a forbidden reaction in SR into a viable reac- 
tion [12, ■ In the case of neutrinos one has the pos- 
sibility of neutrino decay according to reactions such as 
v — > vvD (neutrino splitting), v — > i'e~^e~ (pair produc- 
tion), or — > (photon emission). All these reactions 
are kinematically forbidden in SR, as one can easily see 
by examining them in the neutrino rest frame. However, 
they may be allowed if a neutrino with a given momen- 
tum has an energy greater than the SR energy as it is 
the case for superluminal neutrinos. The corresponding 
decay rate may be written as a function of the neutrino 
velocity, which allows to analyze the compatibility be- 
tween the results obtained from OPERA and SN1987A 
for the neutrino velocity and the experimental observa- 
tion (or non-observation) of the long distance propaga- 
tion of neutrino fiuxes of the corresponding energy. 

In this work we will consider in detail the neutrino 
splitting, since it is a process which only involves the 
neutrino, and one does not need to make further as- 
sumptions on how the departures from SR affect other 
particles. Nevertheless, since the energies that we will 
consider are much larger than the electron mass, the elec- 
tron pair production process will have a similar (order of 
magnitude) rate as the neutrino splitting process and the 
consequences and constraints that both type of reactions 
imply on the difference of the neutrino velocity with re- 
spect to the speed of light will also be similar. Finally, 
the photon emission process, being described by a loop 
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diagram, is suppressed with respect to the other two and 
will give less severe constraints. 



PHENOMENOLOGICAL CONSISTENCY OF 
UNSTABLE SUPERLUMINAL NEUTRINOS IN 
OPERA AND SN1987A 



NEUTRINO SPLITTING WIDTH 

The splitting rate of a neutrino with a Planck scale 
modified dispersion relation 
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at tree level has been estimated in Ref. [12|, yielding a 
decay width 
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where the neutrino energy i?^ is well above the threshold 
of the process, so that the masses can be neglected in the 
calculation. 

The result contains a factor of 3 due to the three pos- 
sible combinations of vv in the final state, a factor G\ 
typical of weak processes, factors (27r) characteristic of a 
three body decay, an factor that can be anticipated 
from dimensional analysis and the last factor involving 
the detailed form of the modification of the relativistic 
dispersion relation. In fact the computation of the neu- 
trino decay width does not depend on the explicit form of 
the modified dispersion relation ([T]). In the general case 
one has 
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If one assumes the standard Hamiltonian relation 
v{p) = dE/dp for the velocity then one has 



2 {v{p) - 1) , (4) 



with 
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The neutrino decay width can then be rewritten in 
terms of the velocity of the decaying neutrino as 
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All the discussion from now on will be based on this 
relation between the neutrino decay width and the dif- 
ference of the neutrino velocity with respect to the speed 
of light. 



From the previous result we can estimate the mean free 
path of superluminal neutrinos as 
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Since this expression depends on v{p) instead of v{p), we 
cannot directly apply the OPERA result or the SN1987A 
bound for the neutrino velocity to obtain constraints de- 
rived from neutrino decay. In order to go further one 
needs to specify the difference of the neutrino velocity 
with respect to the speed of light as a function of the 
momentum p, at least up to the highest momenta for 
which we would like to apply Eq. ([7|. 

The simplest possibility is to consider a momentum 
independent velocity of propagation. But in this case 
there is no neutrino splitting since the kinematic discus- 
sion based on energy-momentum conservation reduces to 
the SR analysis. In fact the expression (jH) docs not re- 
produce this result because the proportionality between 
the typical squared transverse momenta of the neutrinos 
in the final state and the difference of the square of the 
energy and the square of the momentum of the decaying 
neutrino (p^ oc {E"^ —p^)), which is the main approxima- 
tion used in the derivation of the result ([3]), is not justi- 
fied. However the stringent limit on w — 1 from SN1987A 
excludes the interpretation of the OPERA result as due 
to a superluminal momentum independent velocity. 

The next simplest possibility to consider is a momen- 
tum dependent power law for the difference of the veloc- 
ity of neutrinos and the speed of light 
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where A is the high-energy scale associated to new 
physics, p ^ A. In this case one has 
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Let us first examine the compatibility between the 
OPERA and the SN1987A results assuming the power- 
law behavior Eq. up to the energies of the CNGS 
neutrino beam used in the OPERA experiment {{Ey) = 
17GeV 0). 

For n = 1, taking into account that the energies of 
the neutrinos from SN1987A are three orders of magni- 
tude smaller than the energy of neutrinos in the OPERA 
analysis, then the bound on the velocity of propagation 
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for neutrinos in OPERA inferred from supernova (SN) 
observations is v — 1 ^ 2 x 10~^, excluding once more 
the interpretation of the OPERA result as due to super- 
luminal velocities. 

In the case of a quadratic momentum dependence 
(n = 2) a velocity excess with respect to the speed of 
light of2.5xlO-'5 at OPERA leads tow-1 - 2.5x10"" 
at supernova energies and a difference of time of flight of 
neutrinos with respect to gamma rays of around 2 min- 
utes. This implies a spread in the time of arrival of neu- 
trinos of the same order, which is only very marginally 
compatible with the detection of neutrinos in a 15 second 
interval [l4|. In any case, we see that a slightly steeper 
decrease in the velocity from OPERA energies to super- 
nova energies would make the OPERA interpretation of 
superluminal neutrinos compatible with SN1987A obser- 
vations. 

Let us now accept a situation in which the OPERA and 
SN1987A results arc compatible and study the implica- 
tions of Eq. ([?)) considering the simple power-law behav- 
ior Eq. ([8]). In the case of supernova neutrinos, taking 
a typical energy of Ey ~ (10-100) MeV and the maxi- 
mum velocity of propagation consistent with SN1987A, 
Eq. ([ro]) gives a decay length for the neutrino much larger 
than the size of the observable universe. However, the 
mean free path of a neutrino with an energy of 20 GeV 
and a velocity (1 + 2.5 x 10~^) as observed by OPERA 
is ~ 5 (n -I- 1)^ x 10^ km, which implies the impossi- 
bility of the observation of cosmological neutrinos of this 
energy. 

In fact we can use the relation ([T0| to put limits on 
the velocity from different observations. We have 



1 



1 



lOOMeVy 



(11) 



and then lower bounds on the mean free path inferred 
from the detection of neutrinos of a given energy from 
a distant source can be translated into upper bounds on 
the velocity of superluminal neutrinos. In the case of 
supernova neutrinos, for an energy of i? = 20 MeV, the 
bound If, > 50kpc gives (u — 1) < 6 (n-f 1) x 10~^, which 
is less stringent than the one obtained from the time of 
arrival of neutrinos, in agreement with what is said in the 
previous paragraph. In the case of the OPERA experi- 
ment, the bound 1^ > 730 km (the distance between the 
emission at CERN and the detection of the neutrinos at 
Gran Sasso), gives {v—l) < 7 (n + 1) x 10^^, so that the 
simple observation of the expected flux at OPERA puts 
an upper bound on the superluminal velocities that can 
be measured in such an experiment which is only slightly 
above the reported value. 

In conclusion, the OPERA result for a superluminal 
velocity of neutrinos at energies around 20 GeV and the 
bound obtained from SN1987A at an energy three or- 
ders of magnitude lower arc not incompatible with the 



constraints coming from the process of neutrino split- 
ting. However, if one takes seriously the OPERA result, 
these constraints will have deep implications for neutrino 
physics at a higher energy, as we will now examine. 



CONSEQUENCES OF SUPERLUMINAL 
VELOCITIES IN NEUTRINO OSCILLATIONS 
AND COSMOLOGICAL NEUTRINOS 

The decay of superluminal neutrinos may have impor- 
tant consequences in neutrino oscillations and the obser- 
vation of cosmological neutrinos if one accepts the valid- 
ity of the OPERA measurement. From the results ob- 
tained above, it is clear that neutrino splitting does not 
play any role in solar neutrino experiments, which involve 
energies of the same order as for supernova neutrinos. In 
order to consider the implications for atmospheric neu- 
trino oscillations, we need again to make some assump- 
tion on the energy behavior of u(p) at energies around 
and higher than the lOGcV scale. 

In the previous Section we saw how the OPERA re- 
sult can be made consistent with the observations from 
SN1987A if one assumes a dependence of the velocity of 
neutrinos which grows quadratically or quicker with the 
neutrino momentum between both energy ranges. How- 
ever, OPERA also reports an independence of the ve- 
locity in the energy domain explored by the experiment, 
within the statistical accuracy of the measurement 0. 
The consistency of this result with SN1987A makes there- 
fore necessary a change in the momentum dependence at 
this energy domain. 

As an example, let us consider a simple model for v{p) 
satisfying both the compatibility between OPERA and 
SN1987A, and the independence in energy at the OPERA 
energy domain. We consider a "low-momentum" depen- 
dence for the velocity 



v{p) - 1 = (wo - 1) 



(12) 



with n< an exponent which has to be slightly bigger than 
2 covering the momentum interval between OPERA and 
SN1987A. This is combined with a "high-momentum" 
dependence 
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which induces a suppression of the deviations from SR at 
high energies implemented by an exponent n> trying to 
escape to contradictions with high energy neutrino ob- 
servations. The scale Aq which parametrizes the turning 
point from an increase to a decrease of the neutrino ve- 
locity should be in the range of the OPERA recent obser- 
vations in order to reproduce the values for the velocities 
in the high and low energy bins 0. 
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With the previous momentum dependence one can cal- 
culate 



v{p) - 1 = (wo - 1) 
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and the mean free path of "high-momentum" neutrinos 
which results by inserting in Eq. ([7]) the result of v{p) — 1 
in this simple model. In the case p 3> Aq one has 
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Eq. gives then a decay length of the order of 

the Earth diameter for neutrino energies of the order of 
100 GeV. One should expect then a clear suppression in 
the high energy tail of the downside atmospheric neutrino 
spectrum with respect to the upside neutrino spectrum 
as a consequence of the splitting of superluminal neutri- 
nos. Alternatively, the absence of such suppression would 
cast serious doubts on the validity of the OPERA result. 

The consequences for cosmological neutrinos depend 
once more on the evolution of v{p) well above the energies 
explored by OPERA, but in general Eq. ([T6| also implies 
the absence of "high-momentum" neutrinos coming from 
galactic and extragalactic sources. The nonobservation 
of a diffuse neutrino flux on top of the high energy atmo- 
spheric neutrino flux by IceCube [l5| could be due to the 
cutoff induced by neutrino splitting. On the other hand 
a future observation of extragalactic high energy neutri- 
nos will also give very strong constraints on superluminal 
velocities for neutrinos. 



CONCLUDING REMARKS 

A consequence of superluminal neutrinos is the pos- 
sibility of neutrino decay. We have considered the phe- 
nomenological constraints of this fact in relation with the 
recent claim of observation of superluminal neutrinos by 
OPERA. The OPERA resuh and the SN1987A bound 
can be made compatible in the context of departures of 
SR of the form 



with 
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With n< = 2, Ao = 20 GeV and wo - 1 = 2.5 x 10^^ 
one has a scale A for the deviations from SR of 4 TeV, 
many orders of magnitude below the Plank scale which 



is the natural scale for Lorentz violating effects induced 
by quantum gravity. It is this aspect of the recent claim 
of OPERA which makes so difficult to find a framework 
consistent with all observations of neutrinos. The simple 
model considered in this work is a purely phenomenolog- 
ical attempt in this direction. 

We have estimated the decay length of neutrinos for 
neutrino splitting as a function of their velocity. This 
allows us to study the consistency of unstable neutrinos 
with the OPERA and SN1987A results. While both ob- 
servations can be made consistent with the possibility 
of neutrino decay, we find that a value just one order of 
magnitude larger for w — 1 of superluminal neutrinos than 
the value reported by OPERA would make them decay 
before arriving to their detector. 

In our study of neutrino splitting we have neglected 
neutrino masses and the corresponding threshold for the 
reaction. Such a threshold will be pth ^ ^/m^A which 
is of the order of MeV, much below the energies consid- 
ered in the neutrino splitting analysis. This justifies the 
approximation to neglect neutrino masses. 

The decay of superluminal neutrinos may also be rel- 
evant in the observation of atmospheric neutrinos and 
cosmological neutrinos. We have shown that the high- 
energy tail of the atmospheric neutrino spectrum is af- 
fected by superluminal neutrino splitting. The OPERA 
result also implies quite generally the impossibility to ob- 
serve cosmological neutrinos. 

Our analysis rely on the validity of Eq. ([6]) for the 
decay width in neutrino splitting. A limitation of that 
formula is that the approximations made to derive it are 
not valid in the case of a momentum independent veloc- 
ity. This does not affect the analysis made for energies 
lower than the OPERA range, since the compatibility of 
OPERA with SN1987A imposes an energy dependent ve- 
locity for which Eq. ([6|) is correct. However, the OPERA 
data suggest a constant value of the velocity at this en- 
ergy range. If this were a true effect going on for larger 
energies, Eq. ^ would loose its validity and the analysis 
for high energies would change. 

While we were completing this manuscript two new 
papers appeared pointing out the possibility to put con- 
straints jon superluminal neutrinos from neutrino de- 
Our conclusions differ from those of Ref. [idl 
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cay 

due to the different formulae for the decay width. While 
we know that our expression requires a momentum de- 
pendence of the velocity justifying the approximations on 
the angular dependence of the decay, it looses its validity 
for a near independent velocity. In Ref. [l3| the rele- 
vance of departures from SR in the reactions responsible 
for the production of neutrinos is proposed as a comple- 
mentary source of constraints together with the decay of 
neutrinos. A generalization of these analysis including 
a general dependence on momenta for the superluminal 
velocity of neutrinos should be considered. 

As a final word of caution, the present analysis (and 
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those of Refs. [16|, |l7|) rely on the use of the standard 
energy-momentum conservation law, which converts a 
forbidden reaction in SR into an allowed one for superlu- 
minal neutrinos. There are however theoretical scenarios, 
such as Doubly Special Relativity (DSR) theories (isj . 
which consider departures from SR while maintaining a 
relativity principle, requiring then a modification of the 
conservation laws. Quite generally, the consistency be- 
tween the conservation laws of the energy momentum, 
the dispersion relation, and the relativity principle pro- 
duces in these frameworks cancellations which tend to 
hold the SR results of the kinematical analysis of a pro- 
cess 1^, 2^- This could even forbid the decay of supcr- 



luminal neutrinos, depending on the specific DSR model. 
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